Dynamic density functional study of driven colloidal particles: the effect of the system 
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The dynamical properties of classical fluids at pico-liter scale attract experimentally and theoretically much 
attention in the soft-matter and biophysics communities, due to the appearance of the microfluidics ref.[l|, also 
called 'lab-on-a- chip' ref.0, as one of the most relevance 21th century technologies. Theses devices try toput a 
full chemical laboratory onto your hand thought a sophisticated dispositive that includes particles pump ref.Q and 
gates, chemical reaction ref.Q, etc.. In this work we focus our attention on one of the key factors for the design and 
construction of this type of devices: the system dimension. We consider a generic system formed by a dilute solution 
of colloidal particles dragged ( at a constant rate c in the z direction) by a moving potential barrier, modeled by a 
time dependent external potential ( V ex t(r,t) = V ex t(r — ctz)) acting on the colloidal particles but with no effect on 
the solvent. We base our results on a new technique named Dynamic density functional (DDF) theory ref.0,0 which 
is a generalization of Density functional (DF) theory to out of equilibrium, for systems with relaxative molecule 
dynamics (those in which hydrodynamic modes are not relevant), under the hypothesis that the dynamic correlations 
can be approximated by those in a system with the same density distribution at equilibrium; the DDF is a promising 
way of describe classical fluids dynamics at a microscopic level including our knowledge of the very well developed 
DF theory for the equilibrium. 



The central deterministic DDF equation is, 
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where T is the Helmholtz free energy in DF theory, usually split onto the ideal gas contribution and the excess due 
to the interactions (AT), and r o is the mobility of an isolated colloidal particle in the solvent. 

In the case of the shifting external potential we are dealing with, T4 x t( r — ctz) = V e xt( r ') (where r' is the coordinate in 
the reference framework of the external potential), the most relevant feature comes from the stationary state (achieved 
after a transient period from any initial state ) in which p(r, t) = p(r — ctz) = p(r') i.e. when the density distribution 
is shifted at the same rate c as the potential barrier, so that, eq.Q) reduces to, 
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For non interacting particles, such that AT = 0, eq.J2J) becomes a linear Fokker-Planck equation V 2 /0+V(p-V/314) = 0, 
with a "kinetic potential" given by /3T4(r') = f3V cxt (r') + cz' , as a function of the reduced shifting rate c = [3c/T a , 
with inverse length units and (3 — j^p- In this case, the effective dimension behaviour is completely determined 
by the difference between the system dimension (container) and the dimensions over which the external potential is 
homogeneous, so that an effective ID system got thought a 2D plane barrier dragging particles in a 3D container, is 
equivalent to a point-like barrier dragging particles in a narrow channel ref.Q; similarly, an effective 2D system may 
be observed for a cylindrical barrier in a 31? container, or with a point-like barrier in a 2D container. For AT =/= 
(interacting particles) ref . d, this is not longer true, because the excess free energy depends explicitly on the 

dimension of the colloidal particles container. 

Since the system dynamic behaviour is determined mainly through the particles interaction and the effective 
system dimension, and we are interested here in the characteristics induce by this last factor, we consider in this 
work only non interacting particles which would represent diluted solutions of colloidal particles. We set the notation 
before we start solving eq. J2J for each dimension. In the following, we refer to the front (wake) region as that away 
from the external potential, well beyond (behind) the z' positions with V ex t ^ 0; if we assume that this external 
potential barrier is located around the origin of the shifting reference framework, and restricted to values \z'\ w a 
with a the molecular diameter, then the points with z' 3> o are in the front while those with z' <C — o are in the 
wake. In both regions V ext = 0, eq.® reduces to V 2 p + cjp- = 0, and the analytic form of p(r') may be obtained as 
a series with the appropriate asymptotic limit in each dimension. Notice that only for effective ID systems, the front 
and the wake are fully separated by the external potential. In higher dimensions, they are connected through those 
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FIG. 1: Steady state contour density field of non interacting colloidal particles a) for an effective 3D system and b) for an 

effective 2D system, dragging at constant velocity crc = 10 by an external potential V ex t = V e~' r / CT ' with f3V = 10. Dark 
regions correspond to low densities while the bright one show high densities. 



regions with \z'\ < a but |r'| ^> a, for which V ex t — 0. Our result are presented as follow: first we solve analytically 
eq.@ (where V ex t — 0), them we find numerically the full solution, and finally we compare both. 

For a shifting repulsive barrier in effective ID systems rcf. 8], the solution of eq.J2J) is given by p(z') = po + 
Aexp(-cz') at the front, where A/ po is a constant, fixed by the nature of the barrier and the boundary conditions, 
and by p(z') = po, without any kind of structure, in the wake region. It may be proved by a formal reduction of 
the stationary DDF equation with AJ 7 — in any dimension, that this exponential decay behaviour at the front, 
and this absence of excess molecules in the wake, are generic characteristics of the total excess density distribution, 
under constant bulk boundary conditions (lini| r / 1 p = p ), for the integrated excess density distribution across the 

transverse directions, i.e. An ^ z ^ = ~ j dx{^-^- — 1) in effective 2D systems, and AN ^ Z ) — 2j J jj r R(£tLl _ i) m 
effective 31? systems. 

In effective 2D systems, for a shifting cylindric external potential (e.g. the dielectric potential created by a laser beam 
through a bulk solution), eq.J2J), for the region with I4 x t(r') = 0, reduces to, 
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Through a Fourier transform, the solution is given by, 
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p(x, z') — po + I da cos [ai] e @ z g(a); 
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where, 



(5) 



Far away from the external potential, a fast convergence of the Fourier components g(a) is observed, and the 
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relevant features come from their behavior for a <C c. Hence we may use f3+ ~ c + ^- and f3- ~ — for the front 
and wake regions respectively. The expansion of g±(a) as an even polynomial function for small a, and the zero wake 
requirement, lead in the wake to g- (a) « Aia 2 +A2a 4: + ..., while at the front we expect g+{a) ~ Bq + Bio 2 - 
Thus, 
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FIG. 2: Steady state density profile of colloidal particles under the influence of an external potential V ex t = V e~^ r //<T ' , with 
j3V = 10, for effective 2D system at drift velocity ac = 10 (solid line) and ac — 2 (dotted-dashed line) and for effective 3D 
system at drift velocity ac = 10 (dashed line) and ac — 2 (dotted line). The fixed radial and perpendicular distance to the 
z'axis are a) R — x — 0, b) R — x — 0.75a, c) R = x = a, d) R = x = 2a. 



for z' <C — cr; where w = y/c/ \ z' \ and a — wx. 

In effective 3D system, for a shifting spherical external potential (e.g. a driven colloid particle perturbing a bulk 
solution of others colloids), we use cylindrical coordinates r' = (R, <f), z'), and eq.Q, for the region with Kxt(r') = 0, 
reduces to, 



RdR\ BR J dz' 2 dz' 
Through a Hankel transform, the solution is given by , 

/>oo 
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(9) 



where Jo is the zeroth order Bessel function. Proceeding with the Hankel components as in the effective 2D with the 
Fourier ones, we might take f-{a) ~ C\o? + C2C* 4 + ... in the wake, and at the front f+(a) ~ Do + D\o? + D2Q. 4 + .... 
Thus, 



p(R, z) S3 po + e 5z e ~ 
for z' 3> cr, and 

p(i?,z') w + 
for 2' <C — cr; where h — wR. 
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The coefficients (A,B,C and D) are determined by the particular choice of the external potential but, for each 
dimension, the asymptotic decay forms are generic. In the transverse plane (for fixed z') the excess stationary density 
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FIG. 3: Integrated excess of the stationary density distribution over the transverse directions, a) for an effective 2D and b) 3-D 

system. The external potential V ex t = V e~' T , with (3V = 10 is shifting at the velocities ac = 10 ( solid lines), ac = 5 ( 
dotted lines) and ac — 2 ( dashed lines) . 

distribution ( p(r') — p a ) is given by a gaussian times a polynomial function, and along the direction of the movement 
(z' axis) for a fixed distance to the origin (a — const in effective 2D systems and h — const in the 3D ) is given, at 
the front, by an exponential decay (exp(— cz 1 )) and an inverse power in the wake. Over the transverse directions, at 
the front, the maximum density is at the z' axis, while in the wake, this axis is a minimum between two symmetric 
maxima which produce a cap-like structure with parabolic shape. The distance between the maximum and the z' axis 
is 2 y/2 | z' | /ca in effective 31? systems and a/6 | z' \ /ca in the 2D. In fig.l we present the numerical solution for 
the stationary density distribution of colloidal particles under the influence of an external potential V ex t — V r e~' r , 
with /3V = 10, shifting at the velocity ca = 10, in a) an effective 3D and b) 2D system. The clear cap-like structure, 
extended further away from the external potential, is clearly observed. 

We have mentioned before that only in effective ID systems the front and wake are fully separated by the external 
potential, which implicate that as the solution of eq.©, the front and wake should share an unique analytic expression 
in effective dimensions beyond ID. This unique expression produces, at the front (wake), an exponential growth of 
the /?+(/?_) contribution, which only by the appropriate behaviour of g(a)(f(a)) for large a be canceled. We have 
compared our analytical prediction with the numerical solution of eq.J^J), the concordance reveal a local asymptotic 
convergence which makes useful those asymptotic expressions fecis.<|rj)l.(|7l). illOII .( |ll|) ). to represent separately (at the 
front and in the wake regions) the stationary state. The use of polar (spherical) coordinates, instead of the rectangular 
(cylindric) ones, would not solve the problem, because the parabolic structure of the wake implies the entanglement 
of the radial and the angular coordinates. 

Because this work may be usefully in the construction of particles pumps, we want to point out some effects to be 
considered in the design of an efficient device to drag particles. For the same radial form of the external potential 
it is clear that an effective 2D system is more efficient dragging particles than a 3D, because it produces a higher- 
density at the front peak. In fig. 2 we present the numerical stationary density distribution of an effective 3D (dashed 
and dotted lines) and 2D system (solid and dotted-dashed lines), for different fixed values of x and R. The higher 
front density peak in the effective 2D system is clearly observed. For the same reason, an effective ID system is more 
efficient than the 2D. 

For an effective ID system we obtained in our previous work ref.j^] a non-linear dependence of the total current 
(j = c J (p(r') — p )dr') with the drift velocity. If we denoted by c max , the velocity at which the total current is 
maximum for a given external potential, the efficiency of a device as a particles pump increases with c until c max 
is reached, and then diminishes for higher c. As we have mentioned before, the integrated excess of the stationary 
density distribution across the transverse directions, An in effective 2D systems and AiV in the 3D, satisfies the 
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same stationary DDF equation as the ID density distribution, so we may expect for the total current this non-linear 
behaviour too. To support numerically our analytical result in fig. 3 we present the integrated excess of the stationary 
density distribution, a) in an effective 2D and b) 3D system, the exponential decay behaviour at the front and the 
absence of structure in the wake as for the ID stationary density distribution is clearly observed. We have estimated 
that the maximum total currents are obtain for o~c max ~ 0.5 in ID, ps 12.41 in 2D, and « 14.76 in 3D, which 
implicate that, apart from the dependence of c ma x with the external potential shape, it is a function of the effective 
system dimension. c max determines when for a given external potential the system is saturated by the drift velocity. 
Increasing the dimension, we increase the paths for the molecules to avoid the external potential so a higher velocity 
is needed to make the molecules go thought the external potential. 

We have considered the geometrical effects of the effective dimension in systems of non-interacting colloidal 
particles driven by external barriers with constant drift rates. The main difference between the 2D and 3D cases 
analyzed here and the previously reported ID case is the appearance of long wake structures behind the moving 
barrier. In contrast with the common exponential decay p(j') — p ~ exp(-cz') for the structure at barrier front, 
the structure of the wake decays only as an inverse power, p(r') — p ~ l/\z'\ 2 for 3D and p(r') — p ~ l/|z'| 3 / 2 
for 2D. However, the transverse integral of the wake, over the directions other than z' , gives no excess of particles, 
with respect to the unperturbed density p a far from the potential barrier. This is in agreement with the complete 
lack of wake structure behind a ID barrier, and it implies that the global current induced by the shifting potentials, 
calculated as the integral of c(p(r') — p ) comes from the balance between the positive contribution from the 
density peak at the front, and the void (/j(r') < p ) within the barrier. This total current follows a non-monotonic 
dependence with the shifting rates, as reported in ID p|, with an optimum shifting rate c max which increases with 
the effective dimensions. The colloid-colloid interactions effects, which have being leaved out of this work, may be in- 
cluded within the general framework of the DDF 0,0,0] and may help to the design of efficient particle pump devices. 
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